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Abstract.  The  paper  addresses  the  reliability  of  computational  Euialysis  of 
plasticity  problems  in  light  of  the  experimental  work  and  its  mathematical 
formulation.  The  paper  also  addresses  some  basic  constitutive  laws  used  in 
engineering  and  proposes  one  which  leads  to  a  well  posed  mathematical  problem 
and  agrees  well  with  the  experimental  data.  The  data  of  the  experimental 
project  for  5454  aluminum  alloy  in  the  H32  condition  for  periodic  suid  random 
strains  au'e  presented.  Vairious  constitutive  laws  aure  assessed  in  light  of 
the  experimental  data.  Outstanding  problems  related  to  the  reliability  of 
the  analysis  are  mentioned. 

1)  Introduction 

The  problem  of  the  reliability  of  the  computational  analysis  of  solid 
mechanics  relies  on  the  assumptions  made.  In  this  paper  we  will  discuss  the 
question  of  nonlinear  constitutive  laws  especially  related  to  the  plasticity 
problem. 

The  study  of  plasticity  and  cyclic  loading  has  generated  a  large 
amount  of  literature.  We  refer  here  for  exeunple  to  [1]  euid  an  extensive  list 
of  references  there.  We  will  be  especially  concerned  about  proposed  models 
describing  plaistlc  strain  hardening  during  general  cyclic  loading  (e.g. ,  of 
random  chauracter).  These  models  are  either  empirical  or  aire  related  to  some 
mechanical  ideas.  As  examples  we  refer  to  [1-7].  The  models  used  in 
practice  are  usually  one  dimensional,  are  tested  in  very  limited  sets  of 
experimental  data  and  are  extended  to  higher  dimensional  problems  by  some 
mechanical  hypotheses  (e.g.,  v.  Mises,  Tresca,  Hencky,  etc.). 

In  this  paper  we  will  address  the  problem  of  the  selection  of  the 
constitutive  law  especially  in  one  dimension.  Also,  the  generalization  to 
higher  dimension  will  be  mentioned  briefly. 
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2)  The  one-dimensional  model  problem 


The  typical  quasistatic  model  problem  (in  one-dimension)  in  solid 
mechainics  is  the  following:  Let  I  *  (0,1],  W  *  (0,oo)  and  f{x,t), 
f(x,0)  =  0  piecewise  analytic  in  t  and  1*2^1^  every  t,  be  given. 

The  problem  is  to  find  u(x,t)  and  0‘(x,t),  (x,t)  e  IxW  such  that 


(2.1) 

(2.2) 


|^(x,t)  =  f(x,t) 


|^(x,t)  *  c(x,t) 


(2.3) 


<r(x,  t)  ®  A(e(x,  t ) ) 


(2.4) 


u(0,t)  =  u(l,t)  =  0 


(2.5) 

where 


u(x,0)  =  0 


^(x,t)  *  (x,t) 


c  (x,t)  *  (x,  t) 


auid  A  is  an  operator  which  maps  any  function  e(t)  (strain)  continuous  on 
(O.oo),  e(0)  =  0  into  the  function  s(t)  (stress)  which  is  continuous  on 

(0,0#)  with  s(0)  *  0  and  s(tQ),  tjj  >  0  depends  only  on  e(t)  for 
0  ^  t  £  t^.  (Hence  we  assume  that  A  in  (2.3)  is  x-independent ) .  Operator 

A  is  the  constitutive  law  operator.  One  of  the  main  problems  is  to  have 
an  operator  which  guarantees  reasonable  properties  of  the  solution  of 
(2.  l)-(2.5)  such  as  the  existence  and  uniqueness.  In  addition  the  operator 
hais  to  approximate  the  material  behavior  well. 

Let  us  now  formulate  a  specific  class  of  constitutive  operators  A 
which  describe  the  plastic  behavior  and  guarzmtee  that  the  problem 
(2.1)-(2.5)  is  a  well  posed  problem.  We  will  introduce  the  family  of 
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the  constitutive  laws  which  are  based  on  two  following  beisic  assumptions: 

1)  Existence  of  a  convex  yield  surface. 

2)  The  normality  condition:  the  plastic  Increment  is  proportional  to  the 
outwaird  normal  to  the  yield  surface  during  plastic  flow. 


Let 

U  be  a  convex  set  in 

and  let  F(s,a)  be  defined  on 

R  X  11 

so  that 

(2.6) 

F(s,a)  is  convex  and 

functional  on  Rxil. 

(2.7) 

F(0,0)  =  0. 

(2.8) 

There  exist  constants 

jr,  r  such  that  0  <  y  i 

dF 

ds 

» 

dF 

da 

<  r 

uniformly  on  the  set 

{(x,a) |F(s,a)  =  Zq>  for  some 

)  > 

0. 

The  function  F(s,a]  =  is  the  yield  surface.  Let  us  now  derive  the 
constitutive  law  using  the  above  two  postulates. 

Let  (s,a)  lie  inside  the  yield  surface,  i.e.,  F(s,a)  <  Z^.  Then  the 
material  will  be  aissumed  elastic 
(2.9a)  s  =  Ee 

(2.gb)  a  =  0. 

In  (2.9a)  E  denotes  the  model  of  elasticity. 

If  (s,a)  lies  on  the  yield  surface,  i.e.,  F(s,a)  =  then  we  will 
split  the  strain  increment  e  into  elastic  aind  plastic  peirts 

(2.10)  e  «  e^  +  e 

8Uld 

(2. 11)  s  »  Ee^  =  E(e-G^) 


Using  the  normality  condition  we  get  the  following  equation 


(2. 11a) 
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where 


f3£f  ,  r?E _ 2L.1 

(^9aJ  ’aa,J 


Hence,  the  constitutive  law  under  consideration  is 


(2. 12a) 
(2. 12b) 


(s.a)  €  & 


(2. 13a) 


(2.13b) 


where 

(2.14a)  S  =  -^(s, a)  |F(s, a)  <  or  F(s,a)  =  emd  S  o| 

(2.14b)  ?»  =  |(s.a)  [F(s,a)  =  Zq  and  >  oj 

Obviously,  operator  A  is  now  defined  by  (2.13)  auid  (2.14)  for  all  e(t) 
which  su'e  continuous  and  piecewise  differentiable.  (We  assume  e(0)  =  a(0) 
s(0)  =  0). 

The  faunily  of  the  constitutive  laws  (2. 12)-(2. 14)  based  on  the  yield 
surface  F(s,a)  satisfying  (2.6),  (2.8)  will  be  called  an  admissible  plasti 
law  or  briefly  an  admissible  law. 


Now  we  have 


Theorem  2.1.  Let  A  be  admissible.  Then  the  problem  (2.1)  has  unique 


solution  with  c(x,t),  a-(x,t)  piecewise  analytic  in  t  for  every  x  and 

2 

L  ( I )  for  every  t . 

Theorem  2.2.  Let  A  be  admissible.  Then  for  every  mesh  and  linear  element 
the  finite  element  solution  of  the  problem  (2.1)  exists  is  unique  and  is 
piecewise  analytic  in  t.  Moreover,  if  we  let  the  mesh  size  tend  to  zero, 
then  the  sequence  of  the  finite  element  solution  will  weadcly  converge  to  the 
solution  of  the  problem  (2.1). 

Theorems  2.1  and  2.2  have  been  proven  in  a  general  form  in  two 
dimensions  (under  weaker  conditions  than  we  mentioned)  in  [8].  See  also  [9]. 


3)  One-dimensional  engineering  constitutive  laws 

In  engineering  the  constitutive  laws  of  plasticity  are  formulated  in 
different  ways,  usually  in  an  incremental  form  using  the  notion  of  reversals. 
Nevertheless  these  laws  can  be  cast  in  a  form  of  differential  equation.  This 
form  will  be  called  the  standard  form. 

Let  us  mention  some  of  these  laws. 

a)  Bilinear  kinematic  hardening  law.  Here  we  have  one  internal 
parameter  (n  =  1)  and 


(3.  la) 

(3.  lb) 


s  =  Ee 


a  =  0 


for 


or 


or 


Is-k  ^a|  <  s- 


. -1 

s  -  k  a  =  s. 


-s  k  =  s, 


aind  s  ^  0 

and  s  k  0 


(3.2a) 

•  • 

s  »  E  e 

r  .  —  1 

s  -  k  a 

-  S-  8Uld 

• 

s 

>  0 

P 

•  • 

for  - 

-1 

0 

• 

(3.2b) 

a  »  E  ke 

P 

-s  +  k  a 

=  Sq  and 

s 

<  0 
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The  operator  A  based  on  (3.1)  and  (3.2)  Is  independent  of  k.  It 

depends  on  the  model  of  elasticity  E,  the  model  of  plastic  flow  and  the 

yield  stress  s^.  Of  covirse  a(t)  depends  on  k.  Usually  k  =  1  is  used. 

Then,  in  a  physical  interpretation,  this  law  keeps  the  shape  and  size  of  the 

yield  surface  fixed,  and  the  surface  only  shifts  by  the  hardening.  The 

pau'ajneter  a  is  the  position  of  the  center  of  the  yield  surface. 

The  problem  is  whether  we  can  cast  (3.1),  (3.2)  in  a  form  (2. 12)-{2. 14) 

fE-E 

with  F(s,a)  satisfying  (2.6)-(2.8).  Selecting  k*  ,  (3.1),  (3.2) 

will  coincide  with  (2. 12)-(2. 14)  with 


(3.3)  F(s,a)  =  [max  F  (s,a),F  (s,a)] 

1  2 

where  [•,•]  means  a  smoothening  operator  in  the  neighborhood  Q  = 

{ (s,a) lFj(s,a)  =  F2(s,a)>.  (This  smoothening  is  only  formal  to  achieve 
continuity  of  F(s,a)).  In  (3.3)  we  have 


(3.4a)  Fj(s,a)  =  s  -  [^E”]  “ 


(3.4b)  FjCb,.)  .  -s  .  [p|-]  «. 

2 

We  mention  that  the  value  k  =  (E-E  )/EE  is  the  only  value  which  allows  to 

P  P 

cast  (3.1),  (3.2)  into  the  form  (2. 12)-(2. 14) . 


b)  Bilinear  isotropic  hardening  law.  Here  the  stsmdard  form  is  as 


follows 


■ 

V 

1 

1 

0) 

(3.5a) 

•  • 

s  =  Ee 

for  - 

or 

.  -1 

s  -  k  a  * 

0 

Srt 

and 

• 

s 

5  0 

(3.5b) 

a  =  0 

or 

. -1 

-s  -  k  a  = 

0 

^0 

and 

• 

s 

>  0 
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(3. 6a)  s  =  E  e 

P  -1  . 

for  s  -  k  a  =  s-  and  s  >  0 

(3.6b)  a  =  E  ke  “ 

P 


(3.7a) 

(3.7b) 


•  • 

s  =  E  e 
P 

a  =  E  ke 
P 


for 


. -1 

-s  -  k  a 


Sq  8uid  s  <  0. 


Once  more,  the  operator  A  here  does  not  depend  on  k,  while  a(t)  does. 

When  k  =  1  then  an  analogous  physical  Interpretation  as  before  holds. 
The  yield  function  keeps  shape  and  its  center  fixed  and  the  hardening 
increases  its  size.  Once  more  the  proper  selection  of  k,  naunely 
fE-Enl 

k  ®  allows  to  cast  this  law  in  the  form  (2.12),  (2.14)  with 


(3.8) 


(3.9a) 


(3.9b) 


We  get  the  operator  A  which  coincides  with  the  one  defined  by  (3.5)-(3.7). 
Functions  F(s,a)  and  F(s,a)  satisfy  the  condition  (2.6)-(2.8). 
c)  Chaboche  law  [21.  This  law  has  become  recently  popular  and  (see 
Sect. 4)  approximates  the  experimental  results  relatively  well.  Although  it 
is  formulated  in  (2]  in  a  different  form,  it  is  possible  to  write  it  in  the 
standard  form. 
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(3.  10a) 

•  • 

s  =  Ee 

(3. 10b) 

II 

o 

(3. 10c) 

M 

O 

(3. lOd) 

h  =  ° 

(3. lOe) 

O 

II 

e  <  e. 


for 


(3. lla) 

(3.  lib) 


(3. 11c) 


i  _  E(c(a-z)  +  b(Q-R)]  ^ 
c(a-z)  +  b{Q-R)  +  E  ® 


•  _  _ Ec(a-%)  • 

^  ~  c(a-:t)  +  b(Q-R)  +  E 

for  ■ 

6  -  Eb(Q-R)  . 

c(a-z)  +  b(Q-R)  +  E  ® 


(3.  lid)  e.  «  e 

n 

(3. lie)  *  e  -  2  g 


(3.  12a) 


•  _  Etc(a+;i:)  +  b(Q-R)I  • 
“  cU+z)  +  b(Q-R)  +  E  ® 


(3. 12b) 


X  = 


_  Ec(a-*-x)  • 

c(a+z)  +  b(Q-R)  +  E  ® 


(3. 12c) 


R  = 


-Eb(Q-R)  • 

c(a+x)  +  b(Q-R)  +  E  ® 


for 


(3. 12d)  4  =  ;  +  2  5 

n  E 

(3. 12e)  *  ® 

with  the  initial  conditions  s(0)  =  ^{0)  =  R(0)  =  0, 


or 

e  ^  0 

e  ^  0 

®h  =  ® 
emd 

e  >  0 


e  =  e^ 
and 
e  <  0 
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(3. 13) 


e^(0)  =  e„  , 


In  [2]  the  law  is  derived  using  certain  physical  heuristic  arguments.  The 
following  coefficients  zire  input  in  the  Chaboche  law. 
a:  kinematic  coefficient 

c:  kinematic  exponents 

Q:  isotropic  coefficient 

b:  isotropic  exponent 

e^:  yield  strain 
E:  elastic  modulus 

In  contrest  to  the  kinematic  and  isotropic  laws  mentioned  earlier,  the 
Chaboche  law  cannot  be  cast  to  the  form  of  an  admissible  law  as  simple  as 
before.  Nevertheless,  am  admissible  law  which  is  close  to  the  Chaboche  law 
cam  be  formulated.  In  [10]  such  law  (with  n  =  2)  wais  proposed  analyzed  in 
detai 1 . 

d)  The  Babuska-Li  law 
We  have 


(3. 14a) 

•  • 

s  =  Ee 

• 

(3.  14b) 

II 

o 

for  •  (s,a, 3)  e  6 

(3. 14c) 

3  =  0 

• 

9 


(3. 15a) 

(3. 15b) 

(3. 15c) 

(3. 16a) 

(3. 16b) 

(3. 16c) 

where 
(3. 17a) 


6  =  -^(s.a.p)  I  F(s,a,/?)  <  or 
Fj(s,a,0)  *  Zq  and  e  5  0  or 
FgCs.a,^)  *  Zq  and  e  i  o| 
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(3. 17b) 

j>  =  . 

+ 

^(s,a,B) 

1  Fj(s,a,j3)  =  Zq 

and  e  >  o| 

(3. 17c) 

9_  =  . 

^(s,a,B) 

1  F2(s,a,B)  =  Zq 

and  e  <  oj 

and  =  s 

0  y 

is  initial  yield  stress. 

Here 

(3. 18a)  Fj(s.a.0)  =  |a^  -  v^a  +  -  VbQ  0  +  s 

(3.18b)  FgCs.a.p)  =  |a^  +  VaEa  +  ^0^-'i/EQ8-s 

and 

(3.19)  F(s.a.0)  =  [max(Fj(s.a,p).  F2(s.a.8)]*. 

We  see  that  this  law  has  two  internal  paj'ameters  (n  =  2)  emd  depends  on 
the  constamts  a,b,c,Q, E,Zq. 

It  Is  easy  to  see  that  F  defined  by  (3.19)  satisfies  (2.6)  and  (2.7). 
The  property  (2.8)  is  more  complicated  and  needs  some  relations  among 
coefficients  of  input  data.  For  the  detailed  einalysis  of  this  law  we  refer 
to  [10].  In  [10]  also  other  laws  are  proposed  and  analyzed  in  detail. 

4)  Experimental  data 

It  is  obvious  that  the  constitutive  law  influences  essentially  any 
computed  data  (for  example,  data  computed  by  the  finite  element  method). 
Hence  it  is  essential  to  analyze  the  reliability  of  the  assumed  material 
behavior  under  realistic  conditions.  By  this,  it  is  meant  the  behavior  of 
commercial  material  (without  laboratory  treatment)  under  various  random 
strain  conditions.  The  experiments  related  to  the  constitutive  law  have  to 
address  the  following: 

a)  Reproducibility  of  the  stress  response  for  the  imputed  strain. 

b)  Selection  and  quamtltative  assessment  of  an  constitutive  law. 
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For  our  experimental  program  we  selected  5454  aluminum  alloy  in  the  H32 
condition.  The  major  alloying  elements  in  the  workhardening  material  are 
magnesium  (2.754),  manganese  (0.854)  and  chromium  (0.1254).  This  particular 
alloy  is  widely  used  in  marine  applications  and  is  generally  available  in  a 
variety  of  structural  shapes  as  sheet,  plate,  tube,  wire,  rod  and  bair.  The 
normal  mechanical  properties  are  tensile  strength  of  40  ksi,  yield  strength  of 
30  ksi  and  elongation  1054.  Samples  of  the  material  were  obtained  from  two 
suppliers. 

For  the  testing  program  three  different  classes  of  strain  were  used:  a 

triangular  (or  constant  magnitude)  strain  functions  emd  two  different  random 

strain  functions  (using  random  number  generator  for  different  means  and 

stauidaird  deviation).  For  these  cases  seven  tests  were  conducted  with  each 

test  having  a  constant  mean  strain  level  ranging  from  -0.654  to  0.  654.  All 

tests  had  a  strain  range  of  1.054  emd  constant  strain  rate  0.001/sec. 

Replicate  tests  were  concluded  for  each  strain,  thus  total  of  84  data  sets 

were  created.  A  typical  data  set  contained  approximately  40,000  data  points. 

<r  <r 

For  two  samples  of  mezisured  data  <r  the  functions  (p^  ,  4'^^ 

n  =  "-ilti)  - 


*1  *  “'I'h’  * 


were  computed  together  with  the  norms 


||V®’||  =  max  1 4)^1 

15iSN  ^ 


max 

ISiSN 


As  an  error  measure  of  the  difference  the  ratio 
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is  used. 

0 

For  the  measure  of  strain  reproducibility  the  ratio  B  was  computed 
(i.e. ,  two  identical  random  strains  were  compau’ed)  auid  B®  <  1.8%  was 
obtained. 

The  constauits  for  the  constitutive  laws  were  computed  for  every  sample 
by  minimization  of  the  (relative)  difference  between  the  experimental  and 
computed  data  leading  to  minimal  error  measure.  Typical  results  are  given  in 
tables  1  2uid  2.  Here  the  results  are  for  28  periodic  and  28  random  strains 
for  10  loading  loops.  In  Table  1  the  error  measure  B*^  is  reported  for  the 
group  of  periodic  and  random  strain.  Here  the  values  of  the  constants  in  the 
constitutive  laws  are  the  averages  of  the  best  fit  constsuits  for  every  class 
separately.  We  report  the  smallest  (best)  amd  lau'gest  (worst)  error  measure 
of  the  difference  between  the  samples  and  the  law  used  as  well  as  the  average 
error  measure.  In  table  2  we  report  auialogous  data  for  the  combined  set  of 
periodic  auid  random  history  (56  samples).  We  use  the  average  constants  of 
the  best  fits  for  all  (56)  saunples  and  the  average  of  the  best  fit  constants 
for  the  random  strain  (28  samples)  only.  Tables  1  and  2  show  clearly  that 
the  periodic  emd  random  strain  lead  to  different  results.  They  also  show  the 
importauice  of  the  selection  of  the  values  of  the  constants  in  the 
constitutive  laws. 

Table  3  shows  the  results  for  the  error  measure  B*^  for  one  sample 
(FC),  the  span  of  10,  20,  50  and  100  loops  for  the  random  strain.  The 
constants  in  the  constitutive  laws  are  the  average  values  of  28  saunples  of 
random  strain  and  10  loops.  We  report  the  error  measure  of  the  difference  in 
stresses  between  the  experiment  (sample  FC)  (mesm  strain  0.006)  and  the 
stresses  computed  from  the  law. 


13 


Table  1 


28  CYCLIC 

PERIODIC  LOAD  HISTORY 

28  RANDOM  LOAD  HISTORY 

BEST 

AVERAGE 

WORST 

BEST 

AVERAGE 

WORST 

ISOTROPIC 

17.5678*/. 

26.  1709*% 

34.8327*% 

25.9049*% 

36.6324*% 

54.3445% 

KINEMATIC 

13.9703*/. 

21.9894*/. 

29.6227*% 

24.9118*% 

30.3313*% 

37.9495% 

CHABOCHE 

10.0580*/. 

13.4342*% 

19.7299% 

9.7581*% 

15.  1567*/. 

21.6382% 

BABUSKA-LI 

8.7208*/. 

13.0558*% 

17.6285*% 

11. 1595*/. 

17.7304*% 

23.6274% 

Table  2 


ALL  28  CYCLIC  PERIODIC  LOAD  AND  28  RANDOM  LOAD  HISTORY  TOGETHER 


MEAN 

OF  56  REAL 

DATAS 

MEAN  OF 

28  RANDOM  LOAD  DATAS 

BEST 

AVERAGE 

WORST 

BEST 

AVERAGE 

WORST 

ISOTROPIC 

20.2992*% 

32.9285*% 

55.9337*% 

20.2992*% 

36.6324*% 

32.2566% 

KINEMATIC 

15.7110*% 

36.6296*% 

55.7103*% 

19.5086*% 

29.0444*% 

37.9495% 

CHABOCHE 

9.  1939*/. 

14.3332*% 

21.9578*% 

9.0194*% 

14.2968*% 

21.6382% 

BABUSKA-LI 

9.5089*% 

19.3730*% 

36.2526*% 

8.9993*% 

16.5079*% 

23.6274% 
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Table  3 


10  LOOPS 

20  LOOPS 

50  LOOPS 

100  LOOPS 

Ft:  -  HH 

18.6020% 

18. 1069% 

17.6902% 

19.0026% 

FC  -  ISOTROPIC 

41.9331% 

45. 4107% 

44. 1721% 

44. 1659% 

FC  -  KINEMATIC 

27.7171% 

33.2344% 

32.8032% 

32.5333% 

FC  -  CHABOCHE 

21.4401% 

21. 1359% 

24. 1460% 

26.6967% 

FC  -  BABUSKA-LI 

21.3331% 

21.7366% 

26.3657% 

24.7862% 

Table  4 


10  LOOPS 

20  LOOPS 

50  LOOPS 

100  LOOPS 

GL  -  ES 

10.2216% 

10.2265% 

10.6647% 

10.6341% 

GL  -  ISOTROPIC 

29. 1335% 

31.9397% 

34.6639% 

36.8122% 

GL  -  KINEMATIC 

24.5334% 

27.4146% 

30.2475% 

32.4730% 

GL  -  CHABOCHE 

12.7273% 

12. 4393% 

11.9824% 

11.5788% 

GL  -  BABUSKA-LI 

12.3832% 

12.0705% 

11.5576% 

11. 1017% 
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The  first  row  reports  the  error  measure  of  the  stresses  for  two  sajnples 
(FC  and  HH)  (for  the  same  raindom  strain).  This  data  characterize  the 
reproducibility  of  the  response  and  are  related  to  the  uncertainty  of  the 
material  properties. 

Table  4  shows  analogous  results  for  periodic  strain  (once  more  the 
constants  in  the  constitutive  laws  were  computed  from  the  random  strain  for 
10  loops  and  28  samples).  The  error  measure  for  the  difference  between 
two  samples  (GS  and  ES)  (meein  strain  is  -0.004)  is  reported  too.  We  see  once 
more  the  difference  of  the  res{x>nses  for  the  random  and  periodic  load. 
Periodic  load  response  is  obviously  more  predictable.  This  cleau’ly  shows  the 
need  to  make  experiments  not  only  for  periodic  strain. 

Figure  1  shows  the  difference  in  stresses  for  the  seunple  FC  and 
different  constitutive  laws  for  the  random  strain.  Figure  2  shows  analogous 
results  for  the  periodic  strain  (for  the  sample  GL).  Figure  3  shows  the 
strain-stress  relation  for  the  sample  FC  amd  random  constitutive  laws. 

We  refer  to  [11]  for  additional  experimental  results,  statistical 
analysis  of  the  obtained  data  and  also  comparison  with  other  laws.  See  [11] 
also  for  details  of  the  experimental  procedures.  Obviously  unreliability 
(uncertainty)  in  the  material  responses  is  large  auid  hence  reliability 
(uncertainty)  in  computed  (finite  element)  data  is  expected  to  be  leu-ge  too. 
We  will  briefly  address  this  problem  in  Section  6. 

5)  Two  dinensional  constitutive  law 

The  problem  is  in  principle  the  same  as  in  one  dimension.  Here  instead 
the  scalar  strain  and  stress  we  have  to  consider  strain  and  stress  tensors. 
Also  the  internal  parameter  space  is  larger. 
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Figure  3. 


Soiti  il»* 


The  stress  strain  relation  for  the 
and  periodic  strain. 


Because  of  complexity  of  the  experimentation,  one  dimensional  laws  are 
usually  used  together  with  additional  hypothesis  such  as  von  Mises,  etc.  In 
the  case  of  isotropic  material  we  need  rotational  symmetry  of  the  law.  The 
generalization  of  Babuska-Ll  law  (see  (lOl  preserving  rotational  symmetry 

leads  to  the  following  function  F(s,a, B)  where  s  is  tensor  with  the 

3  3  1 

components  s=(s,s,s  )€lR,  a*  (a,,  a„,  a,,]  6  R  and  B  e  R  . 

^  X  y  xy  12  3 

F(s.a,B)  =  f(s-i/ac\/5a)  ♦  ^f^(VRa)  +  |  B^  -  P 

2  2  2  2 
where  f(s)=s  +s  -ss  +3s 

X  y  X  y  xy 


R  = 


4/3 

2/3 

0 


2/3  0  ■ 

4/3  0 

0  1/3 


In  the  constitutive  law,  the  strain  e  hais  to  be  in  the  form  c  ,  c  ,  y  . 

jcx  yy  ^y 

For  details,  analysis,  and  proof  we  refer  to  [10]  where  additional  results 
are  presented. 

6)  Reliability  of  the  computational  analyses  of  plasticity  problems 
Reliability  of  the  computational  analyses  depends  on 

a)  Mathematical  modelling  restricting  the  characteristic  features  of  the 
physical  problem,  available  information  and  aims  of  the  emalysis. 

b)  Correctness  of  the  mathematical  formulation  guarauiteeing  the  basic 
properties  as  existence  of  the  solution,  etc. 

c)  Numerical  treatment  which  produces  aui  approximate  solution  which  is 
sufficiently  close  to  the  exact  solution  in  a  sense  associated  with  the  aims 
of  the  analyses. 

As  h2B  been  seen  in  the  previous  section,  an  importaint  characteristic  of 
plasticity  problem  is  the  large  uncertainty  in  available  information. 
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especially  in  the  constitutive  law.  In  2  and  3  dimensional  problems,  the 
uncertainty  described  by  the  error  measure  B  is  of  order  30%-40%  or  more. 
(In  the  one  dimensional  setting  we  have  seen  in  the  best  ceises  uncertainty  of 
order  20-25%. ) 

The  following  problems  arise: 

1)  How  much  the  of  uncertainty  in  the  constitutive  law  is  projected  in 
the  uncertainty  of  the  data  of  interest? 

2)  Is  it  possible  to  formulate  and  solve  the  problem  to  determine  the 
optimal  bounds  for  the  data  of  interest  depending  on  uncertainties  in 
available  information? 

3)  How  can  one  formulate  and  sol'  '  the  problem  in  a  stochastic  sense? 
Here,  of  course,  various  probability  fields  or  at  leeist  correlations  of  the 
experimental  material  properties  would  be  needed  and  are  practically 
unavailable  today. 

In  [8]  an  attempt  to  derive  bounds  depending  on  the  uncertainties  in  the 
constitutive  law  of  the  worst  scenario  vras  made.  The  idea  is  to  assume  that 
all  uncertainty  in  constitutive  law  is  in  the  unknown  past.  The  problem 
is  then  to  consider  the  "worst"  possible  past  (in  a  set  of  admissible 
function)  leading  to  the  largest  bound  of  data  of  interest. 

There  are  other  possibilities  to  address:  the  influence  of  the 
imcertalnties  auialyzed  by  Monte  Carlo  methods  euid  approaches  of  the  type  used 
in  optimal  design,  etc. 

It  seems  that  any  finite  element  consideration  in  plasticity  problems 
using  constitutive  law  without  addressing  the  effects  of  uncertainties  has 
to  be  interpreted  and  accepted  in  a  very  conservative  way. 

The  problem  of  the  correct  mathematical  formulation  is  a  crucial  one. 

The  formulation  of  the  (direct)  plasticity  problem  guaremteeing  well 
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posedness  can  be  achieved  without  any  adverse  effects  of  approximating 
experimental  data  as  has  been  seen  in  Section  4. 

The  numerical  treatment  should  lead  to  the  approximate  solution  with 
acceptable  accuracy  in  the  sense  needed  for  application.  Today’s  available 
weak  convergence  theory  seems  to  be  far  from  what  is  needed. 

Summarizing  we  see  that  there  are  major  problems  and  questions  related 
to  the  reliability  of  the  computational  analysis  of  the  plEisticity  problems. 
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